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We propose a new orbital controlled model to explain the gate field induced switching of current in a 
semiconducting PbS-nanowirc junction. A single particle scattering formalism in conjunction with a poste- 
riori density functional approach involving hybrid functional is used to study the electronic current; both 
first and higher order Stark effects are explicitly treated in our model. Our calculation reveals that after a 
threshold gate-voltage, orbital mixing produces p-components at the S atoms in the participating orbitals. 
This results in an inter-layer orbital interaction that allows electron to delocalize along the channel axis. 
As a consequence a higher conductance state is found. A similar feature is also found in a PbSe nanowire 
junction, which suggests that this model can be used universally to explain the gate field induced switching 
of current in lead-chalcogenide nanowire junctions. 
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I. INTRODUCTION 

While the size of the conventional silicon based field effect transistor is inching toward its fundamental limit of 
miniaturization, the quantum controlled, semi-conducting nano-wire has emerged as one of the promising options to 
meet the physical challenges imposed by quantum mechanics^i^. Field effect transistors (FET), whose main function is 
to switch source- drain current upon the application of gate field, have been fabricated from semiconducting nano- wires 
(NW) and nano-crystal arrays^"— ; the switching speed for such devices, in some instances, have been found to surpass 
that of the conventional semiconductor-FET. Particularly, IV-VSr— semi conducting nano-wires based FET devices 
have been shown to have huge conductivity and high current gain, which are key requisites for an integrated circuit. 
The controlled synthesis of these NWs with diameter ranging from 1.2 nm to 10 nm have been reported^"— . 

Despite the rapid development on the experimental side, which provides an ample opportunity for theorist to test 
their models, only nominal theoretical efforts are made to understand the quantum phenomenon that dictates current 
modulation in such a nanowire junction (particularly NW of few nm dimension); thus far, no unswerving interpretation 
exists. Since the dimension of the channel is in the nanoscale regime, the electronic property and/or the field effect 
switching mechanism is expected to be different than that of the conventional FET. In this letter, we present a first 
principles quantum transport study in a strongly coupled, single PbS nanowire (PbSNW) junction (Fig. 1) to unravel 
the mechanism responsible for the gate field induced switching of current. We have used the same gold electrode to 
form the nanowire- Zead junction as used in the experiments^— li^. Particularly, we try to answer several fundamental 
questions: How does the gate field affect the intrinsic electronic structure of the nanowire? Can we control the number 
of participating orbitals of the NW- channel via gate field? Which are the orbitals that contribute to the conductance? 
Can we manipulate the shape of the orbital via gating? and finally, is there a universal model that would explain the 
observed gate field induced switching not only in PbSNW but also in other lead-chalcogenide nanowire? 

A single particle scattering formalism in conjunction with a density functional approach is used to study the 
electronic currentii"— . We have included self-consistently both first and higher order Stark effects in our model. Our 
calculation shows, upon application of transverse gate field, the symmetry of the wavefunction is broken along the 
direction perpendicular to the channel axis; the participating molecular orbitals (MO) start to localize in the direction 
of gate field resulting in a shifting of unoccupied energy levels towards the Fermi energy. After a threshold gate voltage 
of -3.7 V, orbital mixing produces a dominant p-component at the S-atom in the participating MOs. This results in 
an inter-layer orbital interaction leading to electron delocalization along the channel axis. As a consequence a higher 
conductance state is found. The higher conductance state is referred here as the ON state and the lower conductance 
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state prior to the threshold value (-3.7 V) is termed as the OFF state. It should be noted that in the OFF state, 
S atom in the participating orbital has only s-component. A similar feature has been observed in a PbSe nanowire 
junction, where s components are found at the Se atom in the OFF state. Thus this orbital controlled model can be 
used universally to understand the observed gate field induced switching behavior in lead-chalcogenide NW junction. 

The rest of the paper is organized as follows. The modeling of the device is described briefly in Section II followed 
by Results and Discussions in Section III. Our main findings are summarized in Section IV. 

II. MODELING THE DEVICE 

For our calculations, we have used a real space approach in which the single determinant wave function is constructed 
from a finite set of Gaussian atomic orbital&i^. This allows us to partition the open NW device structure (Fig. 1) into 
three parts; the first part is the scattering region comprised of a finite NW of length ^ 1.2 nni and diameter ^ 1.17 
nm , the second part is part of the lead that is strongly coupled to the NW and is represented only by a finite number 
of gold atoms (five gold atoms on each side) , and the third part is the unperturbed electrode part which is assumed 
to retain the bulk behavior of gold. The atomic level structural details for the finite part of the PbSNW is taken 
from the optimized structure of an infinite NW, grown in the observed [100] direction. The later structure having 
six Pb and six S atoms in each layer along the growth direction with a lattice parameter of 6 A was calculated using 
the periodic DFT— . Only a five layer NW-structure along the growth direction is considered to build the symmetric 
junction with the lead (Fig. 1). The atomic composition of the lead is taken from the Au [100] surface to avoid the 
lattice mismatch at the NW-lead interface. To realize a strongly coupled junction between the NW and the lead, the 
interfacial distance is varied to determine the optimum distance (2.80 A) where the repulsive interaction is minimum. 

Electron transport is a nonequilibrium (NEB) processii^^^ that requires calculation of the electronic structure of 
the device (Fig.l) under applied bias. The NEB situation refers to the bias condition when the self consistent (SC) 
potential at the lead on one side (Vc) is different from that on the opposite side (V-jz); the equilibrium (EB) situation 
is described by Vc = Vjz. To replicate the NEB situation in our symmetric NW-junction, an electric-dipole interaction 
term is included in the Hamiltonian of the active region (NW+ finite lead) of the device as: Tiled) = 'H(O) +6^ • X)i 
where H(0) is the Hamiltonian in the absence of electric field; is the applied dipole electric field along the axis 
parallel to the direction of current flow (z-axis) , and fl is the coordinate of the i*'' electron; charging effect on the NW 
is considered by including a finite part of the lead. The self-consistent inclusion of dipole interaction term permits us 
to include both first and higher order Stark effects, which is also evident from the comparison of total energy in the 
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active region for different Sd] a non-linear change in energy by increasing the strength of the Sd confirms the inclusion 
of higher order effects. This approach allows us to create an imbalance in charge carrier between the two leads as a 
function of the dipole field strength(ec(); on one lead there is a charge surplus {source) and on the other lead there 
is a charge depletion (drain) resulting local dipoles often referred to as residual resistivity dipoles^^ . This intrinsic 
charge imbalance between the two leads is also reflected from the potential profile summarized in Fig. 2. The relative 
electrostatic potential (REP) in Fig. 2 is calculated by subtracting the average potential at each atomic site in a layer 
at the EB condition from that at the NEB condition. A linear drop in the REP value along the wire axis is noted. 
The magnitude of the potential drop at both the junctions are equal confirming the NW-junction to be symmetric. A 
non linear change in the REP values with different elucidates the nonlinear response of the field. The REP values 
at the left and right Au lead, which are assumed to be at same potential with semi-infinite electrodes on left and right 
respectively, are used to calculate Vc and Vtz- The electro chemical potentials at the two semi-infinite contacts are 
obtained as: fJ^cTZ = ^£,7?. T I-sbT^^^. A small thermal smearing term ksT) in /i£ and /i^ takes into account the 
electronic temperature at the contact in the NEB condition; the potential difference between source and drain [Vsd) 
is then obtained from the difference of /i£ and /it^. In order to simulate the effect of electrostatic gating, we have 
included an additional dipole interaction term (sg • ^he Hamiltonian; the dipole field e'g is applied along 

the direction perpendicular to the channel axis and is referred to as the transverse gate field in this article. In our 
calculation, we have used a posteriori hybrid density functional method (B3LYP) that includes a portion of the exact 
Hartree-Fock exchange. The LANL2DZ effective core potential basis set, which includes scalar relativistic effect, is 
used to describe the Pb and Au atom in the device; a triple zeta augmented by polarization function (6-311G*) basis 
set is used for the S. Subsequently, we recourse to implicit bias- dependent Green's function approacbi2i2i to couple the 
finite NW to the infinite electrode via the finite self-energy functions; coherent, single particle scattering formalism is 
used to calculate the electronic current. 



III. RESULTS AND DISCUSSIONS 
A. Current-voltage Characteristics 

The calculated current-voltage (/sd-V'sd) characteristic for PbSNW as a function of Eg is summarized in Fig. 3. The 
strength of Eg is mapped to the gate potential, Vg, by calculating the potential drop self-consistently between the 
terminal atomic layers of the NW along the direction of gate field. For Vg—Q V, a steady increase in current [Isd) is 
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FIG. 1: (Color online) Schematic representation of a PbS nanowire junction; solid arrows show the direction of the applied 
gate field. 
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FIG. 2: (Color online) Electrostatic potential profile of the NW junction in the absence of gate bias for two different Vsd- 



noted with the increase of Vsd- Changing the Vg from V to -3.7 V, though an analogous linear increase in Igd as 
a function of Vgd is observed, the magnitude of the current is found to be higher ('^ 1.5 times at Vgd ^ 0.7 V) for 
Vg~-3.7 V. A further change of Vg from -3.7 V to -5.6 V reveals a considerable increase in Igd- Changing the Vg from 
-5.6 V to -7.6 V, a non-linear feature in Isd is noted. The current at Vsd— 0-74: is found to be 2.32 times higher for 
V^=-5.6 V than that for V"g=-3.7 V; increasing the negative gate potential from -5.6 V to -7.6 V, 1.53 times higher 
current is found at Vsd '~0.7 V. Thus comparing between Vg—0 V and Vg—-7.6 V, 5 times increase in Isd is found 
at Vsd ^0.7 V. To illustrate this behavior, we have plotted Isd as function of Vg in the inset of Fig. 3; a fixed Vsd is 
used. First a slow increase in current (OFF state) is noted till the value of Vg reaches a threshold value (V^*'') of -3.7 
V. After V^=-3.7 V, a steep increase in current (ON state) is observed with the increase of gate potential resulting 
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FIG. 3: (Color online) Current-voltage characteristics with different Vg for PbS nanowire junctions. The insets show Isd-Vg 
plot for different Vsd- 

in a large change in the slope of Igd — Vg. The calculated ON/OFF current ratio value is found to be 6.28 at Vgd of 
0.54 V between Vg=0 V and Vg~-7.6 V. A similar switching feature is also found in a PbSe NW junction (Fig. 4). It 
should be noted that the ON/OFF current ratio of ^ 3.75 between Vg=0 V and Vg=-8 V at Vsd of 0.5 V is observed 
in a recent experiment, where a single PbS nanowire is used as a channel. The magnitude of Isd reported in the 
experiment is in the nA range in contrast to the ~ fiA current observed in our calculations. Several reasons could be 
attributed to the observed differences in Igd- First, in the experiment the channel length and diameter were 10'^ and 
150 nm respectively, where the diffusive transport could be the prevalent mechanism. In contrast, we have considered 
the channel length and diameter to be 1.2 nm and 1.17 nm for practical purposes. Considering an approximate 
exponential decay in current with the length (Z) for the nanowire used in the experiment e~^'; /3-decay constant), 
we would expect the measured current to be of the order of /lA for a few nm channel length, which has also been 
reported experimentally in single PbSe semiconducting nanowire junction^. In addition, we have considered an ideal, 
defect-free nanowire junction. The magnitude of higher current observed in our calculation is also not surprising 
considering the use of static exchange and correlation potential instead of the true dynamical exchange correlation 
corrected potentiaP^"^^. However, the consistent increase of calculated current upon increasing negative gate bias as 
observed in the experiments^ilS reaffirms on the switching phenomenon replicated here. 

To examine, whether the increase in magnitude of current is due to the use of an implicit orbital dependent B3LYP 
functional approach, we have calculated the current in the same device geometry using different exchange-correlation 
functional; same Gaussian basis sets are used for all the calculations. Our results are summarized in Fig. 5. Though all 
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FIG. 4: (Color online) Current-voltage characteristics with different Vg for PbSe nanowire junctions. The insets show Isd-Vg 
plot for different Vad- 




FIG. 5: (Color online) Current-voltage characteristics of PbSNW for different exchange-correlation functionals at Vg—Q V. 

different functionals (SVWN, PW91PW91, PW91LYP)^^ yield similar current-voltage characteristic, the magnitude 
of the current is found to be much higher (^ 3 times at a Vsd of ^^0.7 V) than that obtained with the B3LYP approach. 
An atomic self interaction corrected DFT scheme yielding a lower conductance than the conventional DFT approach 
has been reported in a molecular junctior^. The inclusion of part of the exact exchange from the Hartree-Fock 
formalism in our posteriori B3LYP approach corrects partly the self-interaction error that occurs in the conventional 
density functional method; it represents a substantial improvement in the right direction as evident from the Vsd — Isd 
curve (Fig. 5). 
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FIG. 6: (Color online) Bias dependent transmission function as a function of injection energy for different gate bias at Vsd ~ 
0.6 V. The Fermi energy is set to zero in the energy scale; dotted lines represent the chemical potential window. Notation: LO, 
LI, L2, L3, and L4 refer to LUMO, LUMO+1, LUMO+2, LUMO+3, and LUMO+4. 

B. Bias Dependent Transmission 

To investigate the intriguing features in gate field induced current and to understand the origin of the field effect 
behavior in PbSNW, we have calculated the bias dependent transmission function as a function of injection cnergy(i?) 
for different Vg(Fig. 6). For brevity, we have only considered Vsd '^0.76 V. First, the increase of area under the 
transmission curve within the chemical potential window (CPW) with the increase of negative gate bias confirms the 
observed increase of I^d with Vg (Fig. 3); the non-linear increase in area explains the change of slope in Isd-Vg plot 
presented in the inset of Fig. 3. Analysis of eigenvalues of Hamiltonian for the NW reveals unoccupied levels (shown 
in Fig. 6) contributes to the conduction. Increasing the Vg, the participating unoccupied eigen-channel shifts in the 
direction of Fermi energy. For Vg—0 V, only LO level contributes to the T{E,V) within the CPW. As Vg increases 
more unoccupied levels move into the CPW, resulting in an increase in the density of states within the CPW. To 
quantify the response of the gate field, we have plotted the Stark shift ( — Eq; i-corresponds to different unoccupied 
levels, £g and cq are respectively the orbital energy in the presence and absence of gate field) as a function of Vg 
for different participating unoccupied levels in Fig. 7. A significant Stark shift has been observed. Different levels 
exhibit different shift, particularly at higher Vg. A closer examination indicates a non linear increase of Stark shift 
(SittiEi + ^Yii jlSijEiEj + . . . ) with the increase of Vg. 
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FIG. 7: (Color online) Stark shift as a function of gate bias. Notation: LO, LI, L2, L3, and L4 refer to LUMO, LUMO+1, 
LUMO+2, LUMO+3, and LUMO+4. A fixed Vsd of ~ 0.6 V is used. 
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FIG. 8: (Color online) Schematic representation of orbital controlled mechanism for the PbSNW FET. In the ON state 
(Vg > V^''), orbital mixing produces p-component at the S atoms resulting in an inter-layer orbital interaction along the 
channel (z) axis. The S-atom has only s-component in the OFF state (Vg < Vg""). 



C. Orbital Analysis 



Next, the natural question to ask is: How does the participating orbitals evolve with the gate bias? Does it have 
any correlation with the observed increase in transmission in Fig. 6? To answer these subtle questions, we have 
analyzed participating MO coefficients in the presence and absence of gate field. As expected, for Vg—OV, the MOs 
are symmetric along the direction perpendicular to the wire axis (negative Y axis). Increasing the gate bias to -7.6 
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V, the symmetry of the wavefunction breaks along the negative Y axis; the participating MOs locaUze in the same 
direction resulting the observed Stark shift (Fig. 7). A close inspection of the MO coefficients reveals that S atoms 
in the participating MO have only s-components in the absence of gate bias. Increasing the gate bias beyond the 
threshold value of -3.7 V, p-components develop at the S atoms due to the strong gate field induced orbital mixing 
(Fig. 8). It should be noted that the Pb atom, which has p-component prior to the application of gate field, does not 
exhibit such orbital evolution. The p-components at the S atoms for the higher gate bias mediates inter-layer orbital 
interaction - allowing electron to delocalize along the channel axis (Fig. 8). This explains unambiguously the origin 
of switching in conductance value observed in Fig. 3. It is worth mentioning that very recently orbital gating has 
been observed in molecular junctions^^. 



IV. CONCLUSIONS 



We present a new orbital-control mechanism to explain the gate field induced switching of current in a semicon- 
ducting PbSNW junction. An implicit orbital dependent single particle Green's function approach that employes 
a self-interaction correction scheme is used to calculate the electronic current. A comparative study using different 
exchange correlation functionals shows a quantitative improvement in the magnitude of current for the self-interaction 
corrected scheme over the conventional DFT. Both first and higher order Stark effects are included in our model. 
The consistent increase of calculated current upon increasing negative gate bias as observed in the experiment, and 
the similar orbital evolution in a PbSe nanowire junction upon application of gate field reassure the validity of our 
generalized model, which can also be used to understand switching of current in other lead-chalcogenide NW junc- 
tions. Thus, the present work may serve as a guiding point in designing orbital-controlled nanowire-FET for potential 
applications in new generation electronic circuit. 
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